Dynamics of Ordering of Heisenberg Spins with Torque — Nonconserved Case. I 
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We study the dynamics of ordering of a nonconserved Heisenberg magnet. The dynamics consists 
of two parts — an irreversible dissipation into a heat bath and a reversible precession induced by 
a torque due to the local molecular field. For quenches to zero temperature, we provide convincing 
arguments, both numerically (Langevin simulation) and analytically (approximate closure scheme 
due to Mazenko), that the torque is irrelevant at late times. We subject the Mazenko closure 
scheme to systematic numerical tests. Such an analysis, carried out for the first time on a vector 
order parameter, shows that the closure scheme performs respectably well. For quenches to Tc, we 
show, to Oie^), that the torque is irrelevant at the Wilson-Fisher fixed point. 
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I. INTRODUCTION 

Interacting systems like magnets and binary fluids ex- 
hibit an ordered configuration at low temperatures, con- 
sisting of coexisting, symmetry broken phases. When 
cooled rapidly from the disordered phase at high tem- 
peratures, such systems take a long time to establish or- 
der, primarily because of the slow annealing of the in- 
terfaces (defects) separating the competing domains. At 
late times, the system organizes itself into a self similar 
spatial distribution of domains characterised by a single 
diverging length scale which typically grows algebraically 
in time L{t) ~ t^^'^ . This spatial distribution of domains 
is reflected in the scale dependent behaviour of the cor- 
relation functions. 

In the last few years, a fairly detailed picture of the late 
time behaviour of the correlation functions has emerged 
iQ. The equal-time order parameter correlation func- 
tion C{r,t) = {(t){r,t) ■ (j){0,t)) is a measure of the spa- 
tial distribution of the domains, and at late times is 
found to behave as f{r/L{t)), where L{t) is the dis- 
tance between defects. The autocorrelation function, 
C(0,ti = 0,t2) = (0(0,0) • (/)(0,t2)), is a measure of the 
memory of the initial configurations, and decays at late 
times as L{t2)^^. The independent scaling exponents z 
and A and the scaling function f{x) characterise the dy- 
namical universality classes at the zero temperature fixed 
point 0. 

The scaling results referred to above, have been ob- 
tained in systems coupled to a constant temperature 
heat bath, where the order parameter (/)(r, t) undergoes 
a purely dissipative (irreversible) dynamics. There is no 
dynamics of the order parameter in the absence of this 
coupling. In general however, apart from dissipating into 
a heat bath, the system may have a hamiltonian dynam- 
ics of its own. Such dynamics can be represented by a 
generalized Langevin equation S, 
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where (/)q, is an A^-component order parameter, is a 
coarse-grained free-energy functional, {•, •} is the Poisson 
bracket, and rja is the noise. The first term on the right 
hand side is reversible, while the second is the usual dis- 
sipative force, with /i = 2 or according as whether the 
order parameter is conserved or not. The noise correlator 
{Vc.ir,t)7ji3{r',t')) = 2TkBT6^f3 {~iVY 5{r -r') 5{t-t') 
is proportional to the temperature T. 

As an example let us consider the dynamics of a bi- 
nary fluid given by the generalized Langevin description 
Eq. (|l|). The relative concentration </) of a binary fluid 
is advected by the velocity fleld v = g/p. The gen- 
eralized Langevin equation for (p involves the reversible 
{(j),ga} 5F / 5ga which reduces to the familiar streaming 
term v • V0 [§ . 

One might classify the dynamics according to the al- 
gebraic structure of the Poisson brackets. For instance, 
the Poisson algebra of the components of the order pa- 
rameter could be of the form {4>a, (j^f^} = c, where c is in 
general a complex number. A recently studied example is 
the dynamics of superfluid ordering of a bose gas ^ . The 
dynamics is written in terms of a complex boson anni- 
hilation fleld ip which obeys the Poisson bracket relation 
— i. On the other hand, the Poisson algebra of 
the components of the order parameter could be a Lie 
algebra {(f>a, 4>f3} = Cq^^i/)^, where Ca/Sj are the structure 
constants. A common example is the dynamics of the 
local magnetic moments of a Heisenberg magnet. The 
components of the magnetic moments (pa {a — 1,2,3) 
satisfy the Poisson algebra, 

{Mr, t), (i)p{r', t')} = tan-i ^lir, t) 5{r - /) 5{t - t') , 

(2) 

where e^^-y is the completely antisymmetric tensor in 
three dimensions and 17^ is the Larmour frequency. 

In this paper, which appears in two parts, we discuss 
the phase ordering dynamics of a Heisenberg magnet in 
three dimensions. Part I of the paper is a detailed discus- 
sion of the nonconserved model, while Part II discusses 
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the conserved ordering dynamics. A section wise breakup 
of Part I follows. After a brief account of the dynamics of 
the model (Section II), we discuss the dynamics of phase 
ordering following a quench to T = (Section III). A 
Langevin simulation and an analytical calculation using 
the approximate closure scheme of Mazenko lead to 
the conclusion that the extra reversible term is irrelevant 
at late times. We subject the Mazenko closure to sys- 
tematic numerical tests and show that it is a fairly good 
approximation at late times. This section is of general 
interest, since this is the first detailed numerical 'test' of 
the Mazenko theory for a vector order parameter. In Sec- 
tion IV, we investigate the dynamics following a quench 
to the critical point Tc and show, using a perturbative 
e-expansion, that the reversible term is irrelevant at the 
Wilson-Fisher fixed point. 



II. HEISENBERG MAGNET AND 
PRECESSIONAL DYNAMICS 

The spins 0q (« = 1, 2, 3) in a Heisenberg ferromagnet 
in three dimensions experience a torque from the joint ac- 
tion of the external field (if present) and the local molec- 
ular field. In response the spins precess with a Larmour 
frequency U.^ about the total magnetic field. Coupling to 
various faster degrees of freedom like lattice vibrations or 
electrons, causes a dissipation in the energy and an even- 
tual relaxation towards equilibrium. 

This dynamics follows from the generalized Langevin 
equation Eq. (^, and the Poisson algebra Eq. (^), 



to obtain the equation of motion in dimensionless form. 



dt 
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The free-energy functional F is taken to be of the Landau 
form, 
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The second term in Eq. (|^) is clearly the torque M x H 
where H = SF/S(p is the local molecular field. 

The free-energy functional is rotationally invariant in 
spin space and so the Poisson bracket term conserves 
the total spin. If the dissipation arises from spin-spin 
interactions, then this will conserve the total spin too and 
so n should be taken to be 2. If however the dissipation 
is a consequence of spin-lattice or spin-orbit interactions, 
then the total spin will not be conserved and so /x = 
0. In Part I we will consider the case of nonconserved 
dynamics, leaving the conserved dynamics for Part II. 

Since the noise correlator is proportional to tempera- 
ture, we may drop it in our discussion of zero temperature 
quenches. We then scale space x, time t and the order 
parameter (p as 
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The dimensionless parameter g = {nL/T)yJrJu is the 
ratio of the precession frequency to the relaxation rate. 
To get a feel for the values assumed by g, let us sei^L 
10^ Hz, P ^ 10^ — 10^" Hz, which gives a range of g ^ 
10-3-10. 



III. PHASE ORDERING DYNAMICS : 
QUENCHES TO T = 

Let us now prepare the system initially in the para- 
magnetic phase and quench to zero temperature. We 
study the time evolution of the spin configurations as 
they evolve according to Eq. (||) . We calculate the equal 
time correlator, 

C(r,t) = (0(r,i).<^(r-Fx,t)), (6) 

and the autocorrelator, 

C(0, ti = 0, t2) = (0(r, ii = 0) • 0(r, h)) , (7) 

where the angular brackets are averages over the random 
initial conditions and space. At late times these correla- 
tors should attain their scaling form 



C{v,t)^ f{r/L{t)) 



c(o,ii =o,^2)~L(^2)" 



(8) 



(9) 



The length scale Lit) is a measure of the distance be- 
tween defects and grows with time as L{t) ~ t^/^. We 
compute the scaling function /(x), the growth exponent 
z and the autocorrelation exponent A by (i) simulating 
the Langevin Eq. (^ and (m) Mazenko's closure approx- 
imation. 



A. Langevin Simulation 

We discretize Eq. (0) on a simple cubic lattice (with 
size N ranging from 4(r to 60^) adopting an Euler scheme 
for the derivatives . The space and time intervals have 
been chosen to be Ax = 3 and Ai — 0.01. With this 
choice of parameters, we have checked that the resulting 
coupled map does not lead to any instability. We have 
also checked that the results remain unchanged on slight 
variations of Ax and Af. Unless otherwise specified, all 
calculated quantities are averaged over 30 uncorrelated 
initial configurations taken from a uniform distribution 
with zero mean. Throughout our simulation we have used 
periodic boundary conditions. 

The scaling of the energy density. 
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s=- j dv{{V<j^{v,t)Y) (10) 

can be used to determine the dynamic exponent z. For 
vector order parameters, e ^ L{t)^'^ at late times, where 
L{t) is the length scale beyond which the field of the 
defect is screened by other defects. Figure 1 shows a log- 
log plot of the energy density as a function of t for g = 
0, 0.5, 1, 2 on a QQ^ lattice. The error bars are smaller 
than the size of the symbols. Upto these times [t — 
16000), there is no evidence of finite size effects. The 
slight curvature seen in the data (especially for larger g) is 
due to finite time corrections. The bold line corresponds 
to A/{t + to) where A and to are varied to give the best 
fit to the data. This shows that the data gathered over 
1.5-decades gives a z = 2, independent of g. 

We next calculate C(r, t) at these late times for dif- 
ferent values of g. Figures 2(a)-(b), are scaling plots 
of C(r,t) versus r/t^/"^ (see Eq. (|)) for g = and 
g — \ on a, 40^ lattice (finite size effects manifest at 
t > 12000). The domain size L{t) extracted from 
C{r = L{t),t) = C(0,i)/2, scales as where the 

exponent is again z = 2 (within statistical errors) and 
independent of g. Figure 3 shows that the scaling func- 
tion f{x) is also independent of g, for 3 = 0, 0.5, 1, 2. 
This scaling function is compared (bold line in Fig. 3) 
with the approximate BPT scaling function for 5 = Q , 
f{x) = (3 7/2 7r)[S(2,l/2)]2F(l/2,l/2,5/2;72) where 
7 = exp(— x^/S) and B and F are the Beta and the hy- 
pergeometric functions respectively. 

At larger values of finite size effects become very 
prominent. This can be seen from the form of the late- 
time C(r, t) for g — b (Fig. 4), simulated on a 50'^ lattice 
averaged over 7 initial configurations. The correlation 
function crosses zero at large r, dips through a min- 
imum, and then asymptotically goes to zero (ofcourse 
/ C{r,t) > 0). It is clear from the figure that at these 
times, C{r/L{t)) for g = 5 would be qualitatively differ- 
ent from the scaling function of Fig. 3. However notice 
that the dip decreases with increasing time. This would 
suggest that the dip might disappear at late times 0, 
and that the resulting scaling function would be identi- 
cal to Fig. 3. In the next paragraph wc will argue that 
this dip is a preasymptotic feature and disappears in the 
scaling limit. This will allow us to conclude that the 
scaling function f{x) is indeed independent of g. 

At very late times, the order parameter field has to- 
tally relaxed with respect to defect cores. Preasymptotic 
configurations typically consist of spin wave excitations 
interspersed between slowly moving defects separated by 
a distance L{t) ^ ^, the size of the defect core. Decom- 
posing (f> into a singular (defect) part (fising and a smooth 
(spin wave) part (psm, we calculate the preasymptotic 
correlation function within a perturbative analysis (see 
Appendix A for details). The computed correlation func- 
tion exhibits a dip at r^/f ~ {l + g'^)/g, which disappears 
algebraically in time (see Appendix). The amplitude of 
this dip increases with increasing g. The dip eventually 



goes away with a relaxation time that scales as (t* is the 
time at which the 5 = correlation function first exhibits 
scaling) , 

T{g)^t4l + g'). (11) 

The crossover time T{g) is estimated to be (taking — 
1000 for the 40^ system) — T{g = 1) = 2000, T{g = 
2) = 5000 and r(g = 5) = 26,000. The crossover times 
for 5 > 5 are much larger than the largest time reached 
in our simulation ! Figure 5 is a plot of pre-asymptotic 
C(r, t) — Csing +Csm at a fixed time where Csing is given 
by the BPT form and Csing takes the form derived in 
Appendix A. Two adjustable parameters related to the 
length scale and the amplitude of the spin wave have 
been tuned to obtain excellent fits to the numerical data. 
Based on these arguments we conclude that the scaling 
fimction f(x) is independent of g. 

We now compute the autocorrelation function 
C(0, 0,t) and extract the exponent A (see Eq. (||)). The 
time t ranges from 4000 — 16000 for the 60'^ lattice (aver- 
aged over 50 initial configurations), well into the scal- 
ing regime for the one-time correlator. Figures 6 are 
log-log plots of C(0, 0, t) versus t for various values of 
g = 0, 0.5, 1, 2. It is difficult to give a precise value of 
the decay exponent A, since as can be seen from Fig. 7, 
the "effective" A varies by about 3% over half a decade. 
However it is evident from the bold line in Fig. 6, which 
corresponds to A{t + to)^^^^ {A and to are varied to give 
the best fit to the overall data), that the value of A is 
independent of 5. A fit to each data set for a given 
g, obtains the following values for the exponent A — 
X{g = 0) = 1.526 ± 0.007, X{g = 0.5) = 1.521 ± 0.008, 
X{g = 1) = 1.55 ± 0.01, and A(g = 2) = 1.55 ± 0.02. The 
numerical values listed above can be compared to the 
Mazenko closure estimate of 1.587 for 5 = ||l[. These 
values obey the Huse-Fisher bound A > d/2. 

It is clear that finite size effects set in at later times. 
As discussed in S , finite size effects will be relevant when 
the spread in C(0, 0, t) (given by AC(0, 0, t) - N-^^^) is 
of order C(0,0,t) itself. This will happen when L^^ ^ 
The fact that the numerically computed A in- 
creases marginally with g, indicates that finite size effects 
are more apparent for larger g. This is consistent with 
our discussion on the effects of finite size on single-time 
correlators. 

We end this section with the following assertion based 
on our careful numerics. The exponents z and A and the 
scaling function f{x) are seen to be independent of the 
torque g. This implies that the torque is irrelevant to 
the late-time dynamics at T = 0. In the next section, we 
will apply the approximate method of Mazenko to this 
problem and arrive at the same conclusion. 

B. Application of the Mazenko Closure Scheme 

Of the variety of approximate schemes devised to eval- 
uate the form of the scaling function, the closure scheme 
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introduced by Mazenko [Q is amenable to systematic im- 
provement [yj^. We shall use this closure scheme to de- 
termine the scaling form /(x) and show that it is inde- 
pendent of g. The method consists of trading the order 
parameter </)(r, t) which is singular at defect sites, for an 
everywhere smooth field m(r, t), defined by a nonlinear 
transformation. 



(12) 



At late times, the magnitude of (f) saturates to its equi- 
librium value almost everywhere except near the defect 
cores. This suggests that the appropriate choice for the 
nonlinear function a is an equilibrium defect profile, 
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V2„a(m(r,t)) = r (a(rn(r,t))) 



(13) 



where V'{x) = —x + [x ■ x)x. This choice allows for a 
natural interpretation of m (in the vicinity of a defect) 
as the position vector from a defect core. The simplest 
nontrivial solution of Eq. ( |l3| ) is the hedgehog configura- 
tion, 



CT(m(r, t)) 



m{r, t) 
|TO(r, t)\ 



.9(MI), 



(14) 



where f;(0) = and g{oo) ~ 1. Equation (||) can be 
used to derive an equation for the correlation function 
CY12) = {${vi,ti) ■ ${r2,t2)). Substituting for cj) (Eqs. 
(O) and (M)) in the right hand side of the resulting 
equation, we get 

dtC{l2) - V?C(12) + ( <7(77i(2)) • V'{cr{m{l))) ) 

+ g ( <T(m(2)) • ct(to(1)) x v2<t(to(2)) ) . (15) 

So far no approximation has been made, but further 
progress seems impossible without one. Now along with 
Mazenko, we make the assumption that each component 
of m(r, t) is an independent gaussian field with zero mean 
at all times. This implies that the joint probability distri- 
bution P(12) = P(to(1),to(2)) is a product of separate 
distributions for each component and is given by 1 1| , 



The joint distribution has been written in terms of the 
second moments ^0(1) = (m£j(l)^), Sq{1) = (771^^(1)^) 
and Co(12) = (TO„(l)m„(2)). 

With this assumption, the right hand side of Eq. (|T5| ) 
simplifies to. 



^=V^C(12) + -^^ 
dh ^ ^ 25o(l) 97 

+ g {a{m{2)) ■ a{m{l)) x V'^a{'m{2))\ 



(17) 



It is clear that since P(rn(l), m(2)) = P(— to(1), — m(2)), 
the last term, which is odd in m, drops out. The result- 
ing equation is identical to the purely dissipative one, 
showing that the torque is irrelevant at late times. This 
conclusion, consistent with our earlier numerics, is a di- 
rect consequence of the Mazenko closure approximation. 



C. Justification of tiie Mazenlto Approximation for 
vector order parameters 

In this section we will subject the above stated as- 
sumptions of the Mazenko closure scheme to systematic 
numerical tests. We will do this by numerically solving 
the Langevin equation Eq. (^) by the method outlined in 
Section III A. Knowing (j){r,t), one can compute m{r,t) 
by inverting Eq. (p^. This is facilitated by choosing 



+ m ■ m 



(18) 



which is consistent with the boundary conditions for g{x) 
mentioned in the previous section. The resulting ansatz 
for (T, 



CT(m(r,t)) = 



\/l + rn ■ rri 



(19) 



can be easily inverted. We calculate both the single point 
probability distribution P(m(r, i)) and the joint proba- 
bility distribution P(12) in the scaling regime and com- 
pare with the Mazenko assumption. 



cxp ■ 



1 / ml{l) , m2(2) 27m„(l)TO„(2) 



2(1-7^) \^5o(l) ^o(2) ^Soil)So{2) 

(16) 



where 



and 



27r^(l-72)5o(l)5o(2) 
Co(12) 



7 EE 7(12) 



V^o(l)^o(2)' 
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In what follows all probability distributions have been 
computed on a 45'^ lattice and averaged over 100 initial 
configurations. We shall display the distributions for g = 
and (? — 1. We have collected data in the scaling regime 
from t — 2000 to t = 15000, after which finite size effects 
set in. 

Figure 8 and 9 are scaling plots of P{mi{r, t)) at g = 
and g = 1 respectively. In accordance with the Mazenko 
assumption, the scaling variable has been taken to be 
mi/y%{t), where So{t) = (mi(r,t)2). A plot of So{t) 
(for both g = and 1) versus t (Fig. 10), shows a linear 
growth over a decade, consistent with a z — 2. The scaled 
distribution P{mi) is also seen to be the same for g = 
and 3 = 1 (Fig. 11), suggesting that it is independent of 

g- 

Though the distribution grossly resembles a gaussian 
at late times, a closer inspection shows systematic devi- 
ations at small values of mi (Fig. 11). The distribution 
seems to be flatter than a gaussian when toi « 0. This is 
clearly visible in a plot of — ln( — In [P{m1/ Sq)] ) versus 
ln(m^/S'o), which shows that distribution deviates from 
a gaussian for small mi (Fig. 12). These findings are 
consistent with a similar analysis done on a scalar order 
parameter pX| ]. 

We now study the two point distribu- 
tion P(m(l), m(2)). Equation ( pj[ ) implies that in the 
variables 

the distribution P(12) = P(TO+(12))P(m_(12)), where 



P(m±(12)) = JJVA^exp. 



(l±7)V^o(l)^o(2) 



(20) 



We compute P{ma+{l2),ma-{l2)) a.t ti ^ t2 = t in 
the scaling regime. A plot of P(mi+,mi_) (Fig. 13) for 
|ri — = 4\/3 looks like a product of two gaussian 
distributions. To check this we compute the difference 

A(m„+(12),mo_(12)) = P(m„+(12), m„_(12)) 

-P(m„+(12))P(m„_(12)), 

which should be zero everywhere if the Mazenko approx- 
imation were to hold. 

Figure 14 shows a surface plot of the difference 
A(mi+(12), mi_(12)), magnified 10^ times. It is clear 
that A(toi+(12), mi_(12)) is zero everywhere except in 
the region close to the origin, where the maximum de- 
viation from is around 10^^ . The situation is sim- 
ilar for g ~ 1. We now compute P(mi+(12)) and 
P(mi^(12)), and find that the scaled distributions are 
the same for g = and g — 1. Figures 15 and 17 
are scaling plots of P(mi_|_) and P(toi_) respectively 
and indicate that the scaling function is independent of 
g. Moreover the distribution looks like a gaussian, in 



accordance with the Mazenko theory. Figures 16 and 
18 are plots of - ln(- ln[P(m2^(12)/(m2^(12)))]) versus 
ln(TOf_(_(12)/(mf_|_(12))). The deviation from a straight 
line when mi± « 0, indicates that the distributions dif- 
fer slightly from a gaussian. Note that the data for small 
mi± in Figs. 16, 18 do not quite scale and so it is likely 
that in the computation of the joint-probability distribu- 
tion, we have not yet reached the scaling regime. 

In conclusion, we have shown, by computing the sin- 
gle and the two -point probability distributions, that the 
Mazenko closure scheme is a fairly good approximation 
at late times with possible deviations when nia ~ 0. This 
study justifies the use of the Mazenko approximation 
in the evaluation of the correlation functions and hence 
the results of the previous section. This constitutes the 
first serious numerical check on the assumptions of the 
Mazenko scheme for vector order parameters. 



IV. PHASE ORDERING DYNAMICS : 
QUENCHES TO T = Tc 

We now quench from the high temperature paramag- 
netic phase to the critical point Tc, and ask whether 
the spin precession changes the late time dynamical be- 
haviour. More precisely, we would like to know whether 
the torque term g is relevant at the Wilson-Fisher fixed 
point corresponding to the pure g = Heisenberg model 
with nonconserved dynamics. 

This can be done by power counting to C(e^) at 
the Wilson-Fisher fixed point r* = -(5/22)A^e, u* = 
87r^e/ll (A is the ultraviolet cutoff). Dimensional analy- 
sis provides the scaling dimension [g] ^ d/2 + l~ z + ri/2. 
At the Wilson-Fisher fixed point z = 2 + C77, where 
c = 61n(4/3) - 1, and r? = (5/242)e2. This implies that 
above d = 2, the torque term is irrelevant. In c? = 1, 
the torque is relevant, as can be seen by an explicit so- 
lution of the Langevin equation Eq. (|g) in one spatial 
dimension. 



V. CONCLUSIONS 

We have investigated the effect of reversible Possion 
bracket terms in a generalized Langevin equation on the 
phase ordering dynamics at late times. The dynamics 
of such systems can be classified by the structure of the 
Poisson brackets. In this Part I of two parts, we have 
made a detailed study of the nonconserved dynamics of 
an order parameter whose components obey a Lie Alge- 
bra. A common example is the nonconserved dynam- 
ics of a Heisenberg magnet. The dynamics consists of 
two parts — an irreversible dissipation into a heat bath 
and a reversible precession induced by a torque due to 
the local molecular field. For quenches to zero tempera- 
ture, we have shown, both numerically (Langevin simula- 
tion) and analytically (approximate closure scheme due 
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to Mazenko), that the torque is irrelevant at late times. 
The Mazenko closure was subject to critical numerical 
tests and was shown to perform well at late times (apart 
from small deviations). For quenches to Tc, we show to 
0{t^), that the torque is irrelevant at the Wilson-Fisher 
fixed point. 

In Part II of this series we investigate the effect of 
the reversible torque on the conserved dynamics of the 
Heisenberg magnet in three dimensions. We will show 
that the torque is relevant both for quenches to T = 
and T — Tr- 



APPENDIX A: 

In this Appendix we study the effect of spin wave ex- 
citations on the time dependent preasymptotic equal- 
time correlation function. We would like to show that 
inclusion of such excitations in the correlation function 
C(r, t), leads to a dip at preasymptotic times when g ^ 
which eventually relaxes. As suggested in Section III 
A, at very late times, the order parameter field has to- 
tally relaxed with respect to defect cores. Preasymp- 
totic configurations typically consist of spin wave exci- 
tations interspersed between slowly moving defect cores. 
These defect cores are separated by a typical distance 
L{t) 3> the size of the defect core. In general one 
can decompose (l){r,t) — 0sing(r, i) -I- (/)sm(r, t), where 
the singular part (l)sing parametrizes defect configurations 
while the smooth part (j^sm is a linear combination of spin 
waves of wave- vector k, (/)sm(r, t) = Y^^l"^ ^ (/)^™(i)e**'''". 
The preasymptotic correlation function will thus have 



three contributions 



Cam = {(t>sm,{0,t) ■ (j3sm{'^,t)) and the scattering of spin 
waves from slowly moving defects Cscat = {(f'sm{0,t) ■ 
4>singi^,t)) ■ At late times of course (f>^''{t) 0, and 
(losing can bc traded off for the auxilliary field m within 
the Mazenko approach. Thus the Csmg (r, t) part of the 
correlation function is given by the solution of Eq. ( p^ ) 
or the BPT form §. 

The smooth part of the correlation function Csm (r, t) 
can be estimated from a perturbative calculation, 
wherein the defects separated by a distance L ^ ^ arc 
taken to be static (justified post priori). We shall see 
that the dip is a result of this smooth part. Confining 
our attention to a single domain of size L, we can split 
the smooth into transverse and longitudinal compo- 
nents about the well-defined broken symmetry axis taken 
to be along a = 3. Thus (f>^{r,t) = (peg Sad, + Ua{'r,t), 
where the equilibrium magnetisation (p^q — 1. 

Consider an initial smooth localised pulse in the in- 
terior of this domain of the form ui(r,0) = W2(r, 0) = 

the widths cr, w <C i and u(0), uafO) <C 1. The equation 
for Ua can be read out from Eq. (0), 



~dt 



-{upUp)5a'i " 2(UqU/3)(5/33 
+gea0jUpV'^Uj - UpUfjUa- 



(Al) 



To solve this equatio n pe rturbatively, we multiply the 
nonlinear terms in Eq. (Al) by an arbitrary real param- 



eter e (< 1) and express Ua{v^t) as e^u^^\r,t). The 

initial conditions for ?ii"-*(r, 0) follow from Ua{r,Q). We 
are interested in solutions that decay as t ^ 00. Conver- 
gence of the perturbation series at e = 1 is guaranteed by 
the smallness of the initial deviation and because higher- 
order terms in the expansion decay faster. To ©(e"), the 
spin waves do not interact and 



(0) 



~dt 



V^u<^:-^-gea,^V^uf~2uf5. 



q3 



(A2) 



The equations for u^^^ and u^2^ decouple in the vari- 



(0) 

± 



{uf^ ± iu'^')l2^ giving rise two precessing 



,(0) 



ables u 

goldstone modes in the transverse direction and an ex- 
ponentially decaying mode in the longitudinal direction. 
Thus, 



and 



4°)(r,i) 



(1 ±i)u(0) 
16(7rt(l±i5))3/2 

"3(0) 



exp 



4(1 ± ig)t 



exp 



At 



- 2t 



(A3) 



are the asymptotic solutions to 0{e'^). 

To C(e), the dynamical equations in the transverse 
variables are given by. 



dt 



(1 ±i.g)V2u£^ -2u. 



± ig{u. 



-ufv^u^i^)-{Aufu'-^^ + {uff)uf (A4) 

The last two terms are subdominant in 1/t and 
u(0),M3(0) respectively, and so the transverse correlator 
to 0{e) is given by, 

Csm (!■,*) = ^(w-^(x,<)M_(x + r,i) +u+(x + r,t)u_(x,i)) 



^3/2 



exp 



X < cos 



r 

gr- 



A2e- 



U{3 + .92) 



t3 



-|-7r/4 



■ exp 



2i(3 + .92) 

(A5) 



where Ai - 0(^(0)^) and A2 ^ O(w(0)^ 1*3(0)) are con- 
stants depending on g and initial conditions. The cosine 
term in the above expression results in the observed dip 
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of the total correlation function. The magnitude of the 
dip increases with increasing g. 

The dynamical equation for the longitudinal compo- 
nent to 0{e) is likewise given by, 



[7] 



9u' 



(1) 



dt 



(A6) 



As before, the terms proportional to {u^^^Y, the gra- 
dient terms and the cubic term are subdominant, and so 
the decay of the longitudinal correlation function cjm, is 
given by, 



CL(r,t) = (03(x,i)03(x + r,i)) 



[8] 
[9] 
[10] 

[11] 



^3/2 
't9/2 



■ exp 



exp 



■ exp 



2i(3 



4^(1 + 



.9^) 
(A7) 



where Bi - O{u3{0)^ ,u{0)^), B2 - O{u{0f ^3(0), u{Of) 
and ,63 O(u(0)'*) are funtions of g and initial condi- 
tions. 

Note that C^„ evolves with a width that scales as t^^^ 
and whose amplitude decreases as t~^/^. The longitudi- 
nal Csm decays exponentially fast. This decay is consis- 
tent with our earlier assertion that the defects separated 
by a distance L{t) hardly move over time scales corre- 
sponding to spin wave relaxation. 

The cross correlator Cscat coming from the scattering 
of spin waves by moving defects, can be calculated by 
treating 0sm as "slaved" to (j)sing- As the defects move 
they excite spin waves which decay in a time scale smaller 
than the time taken by the defects to move any further. 
The domin ant contribution to Cscat comes from the prod- 
uct of Eq. ( AS ) and (fusing- It is easy to see that this term 
leads to the same cosine dip as in Eq. (A5) but with an 
amplitude that decays algebraically in time. This is the 
source of the slow decay of the dip. 



In order to make such a claim, great care should be taken 
to avoid finite size effects. The tail of the autocorrela- 
tion function gets afi^ected by the finite size of the lattice, 
earlier than the part of C(r, t) used to calculate z. The 
decrease in the dip at g = 5 can only be perceived when 
the system size is large. 
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1 




0.2 




X = r/L(t) 



Fig. 1 log-log plot of the energy density e versus t for 
various values of g, g — 0(o), g — 0.5(n), g — 1(A), g = 
2(*). The straight line is a fit (see text). 



0.6 




X = r/L(t) 

Fig. 2(a). The scaling function f{x) versus x = r/t^/'^ 
for 5 = 0. Data taken at t = 1000 (o), t = 5000 (A) and 
t — 10000 (★). Error bars are smaller than the size of 
symbols. 



Fig. 2(b). f{x) versus x = r/t^/"^ for g = I- Data 
taken at t = 2000 (o), t = 5000(A) and t = 10000 (*). 
Error bars as in Fig. 2(a). 




X = r/L(t) 

Fig. 3. f{x) versus x = rjt^l'^ for different values of 
5 (g = (o), 0.5 (□), 1 (A), 2 (★)). The data for g^2 
was obtained from a simulation on a 50^ lattice aver- 
aged over 7 configurations. The continuous curve is the 
approximate analytical form defined in the text j^]. 
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Fig. 4. The correlation function C{r,t) versus r for 

g = 5 at various times ( t = 3000 ( ), t = 5000 (- • 

— •), t = 7000 (•••)• Note that the dip gets smaller as 
time progresses. Error bars as in Fig. 2(a). 




t 



Fig. 5. Prcasymptotic C(r, t) for g = 1 and 5 ( bold 
lines ) calculated in Appendix A. Note the excellent fit 
to data for 5 = 1 (o), g = 5 (□) 



Fig. 6. Log-Log plot of the autocorrelation function 
C(0,0,i) versus t for .9 = 0(o), 0.5(0), 1(A), 2 (★). The 
straight hne is a fit A{t+tQ)^^/^ , where A has been chosen 
to be the Mazenko value —1.587. 



-1.51 I < < < > I > < I < I < I < < I < < I < I 




_j 55 I . . . . I . . I . I . I . . I . . I . I . . . . I I 
5x10-5 0.0001 0.00015 0.0002 0.00025 

1/t 

Fig. 7. Effective A as a function of 1/t for g = 0. 
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Fig. 8. Single point distribution of P{m\) for g = 
at different times t = 5000 (o), t = 10000 (□), t = 
15000 (A) . The distribution scales in the variable mi/L, 
where L{t) = ^/S^. 




m,/L(t) 



Fig. 9. Same as Fig. 8, but for g — 1 for times t = 
5000 (o), t = 10000 (□), t = 15000 (A) . 



Fig. 10. Linear growth of So{t) with t for g = (o) 
and g = 1 (•). 




-4 -2 



mi/L(t) 



Fig. 11. Scaling plot of P{mi) is independent of g. 
Data taken for 5 = (i = 5000 (o), t = 10000 (□)) and 
g=l{t = 5000 (A), t = 10000 (*)). 
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1.0 




Fig. 12. Plot of -ln(-ln[P(mf/5o)] ) versus 
ln(m?/5o) for g = {t = 5000 (o), t = 10000 (□)) and 
g=l{t = 5000 {A),t = 10000 (★)). The line with slope 
— 1 drawn for comparison, highlights the deviation of the 
data from a gaussian at smaller mi. 



Fig. 13. Joint probability 

distribution P(to+(12), m+(12)) for 5 = at t = 5000 
and |ri — = 4-\/3. 
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Fig. 14. Surface plot of A(mi+(12),mi_(12)) (magni- 
fied 10^ times) for g = at t = 5000 and |ri - r2| = 4\/3. 



Fig. 16. P(m+i/L+) shows deviation from gaussian 

for small m+. Data shown for g = (t = 5000 (o), t = 
10000 (□)) and £/ = 1 (i = 5000 (A), t = 10000 (★)). 




Fig. 15. Scahng plot of P{m+i/L+), where = 
^{m\^{Y,t)), for g = = 5000 (o), t = 10000 (□)) and 
g = l[t = 5000 (A), t = 10000 (★)). 



a. 




m,_/L_(t) 



Fig. 17. Scaling plot of P(m_i/i_),where L_ = 
sJ{m'i^{Y,t)), for £f = (t = 5000 (o), f = 10000 (□)) 
and g = 1 (t = 5000 (A), t = 10000 (*)). 
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ln(mt/<m;_)) 



Fig. 18. P{m-\/ L-) shows deviation from guassian 
for small m_i. Data shown for g = (i = 5000 (o), t = 
10000 (□)) and = 1 = 5000 (A), t = 10000 (★)). 
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